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Abstract 

We present a high statistics study of the D- and I?-meson decay constants. The 
results were obtained by using the Clover and Wilson lattice actions at two different 
values of the lattice spacing a, corresponding to (3 = 6.0 and 6.2. After a careful analysis 
of the systematic errors present in the extraction of the physical results, by assuming 
quite conservative discretization errors, we find fr> s = 237±16 MeV, fjj = 221±17 MeV 
UdJJd = 1-07(4)), f Bs = 205 ± 35 MeV, f B = 180 ± 32 MeV (f B Jf B = 1.14(8)), in 
good agreement with previous estimates. 



1 Introduction 

Jb is a relevant parameter in phenomenological studies of the Standard Model, in 
determinations of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
and in studies of B—B mixing. 

The leptonic decay width of the B meson is given by 

T(B+ -> r V) = ^^Mb (l " ^1) 2 Mlfl ■ (1) 

A precise knowledge of the leptonic decay constant fs, analogous to f n in tt — ► fiu^ 
decays [], would then allow an accurate extraction of the CKM matrix element \V u b\, 
which is actually known with a relative error of about 25%. 

fs also enters in phenomenological analyses of the B°-B° mixing amplitudes, to- 
gether with the so-called renormalization group invariant i3-parameter Bq. The square 

of the mixing parameter £ = Jb\J Bb is in fact related to the matrix element of the 
renormalized AB = 2 Hamiltonian |ffl|. All theoretical calculations of Bb tend to give 
values very close to one for both the B% and the B® mesons. Thus the strength of the 
mixing is essentially regulated by the meson decay constant. It was realized a long 
time ago j2j that a value of £ ~ 200 MeV, combined with a large value of the top 
mass, leads to a large value of sin 2/3, the parameter which controls CP violation in 
B — ► J/tpK s decays 

N{B« -. Jli>K,){t) + JV(B° -. Jli,K,){t) ' 

Thus a precise theoretical determination of fs is very important. 

One of the most precise determination of fs is presently given by lattice QCD 
calculations. Leptonic decay constants are usually computed from the matrix element 
of the lattice axial current: 



Z A (0 14, | P(p = 0) ) =if P Mp, (3) 

where Mp is the mass of the pseudoscalar meson, 4 = Q(z)jol5Q(x), where Q and 
q denote the heavy and light quark fields respectively, is the fourth component of the 
lattice (local) axial current, and Za is the renormalization constant necessary to relate 
the lattice operator to the continuum one ||, Q. In the case of /„-, fx, and more 
recently fp> a , this is an example of a simple quantity for which a comparison of lattice 
results with experimental data is possible (we notice en passant that fu s was predicted 
by lattice calculations long before its measurement 

The major sources of uncertainty in the determination of fp, besides the effects due 
to the use of the quenched approximation, come from the calculation of the constant Za 

We use the normalization convention in which /„. = 132 MeV. 
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Run 


C60 


C62 


W60 


W62a 


W62b 


P 


6.0 


6.2 


6.0 


6.2 


6.2 


Action 


sw 


SW 


Wil 


Wil 


Wil 


# Confs 


170 


250 


120 


250 


110 


Volume 


18 3 x 64 


24 3 x 64 


18 3 x 64 


24 3 x 64 


24 3 x 64 


Ki 




0.14144 




0.1510 






0.1425 


0.14184 


0.1530 


0.1515 


0.1510 




0.1432 


0.14224 


0.1540 


0.1520 


0.1520 




0.1440 


0.14264 


0.1550 


0.1526 


0.1526 




0.1150 


0.1210 


0.1255 


0.1300 


0.1300 




0.1200 


0.1250 


0.1320 


0.1350 


0.1350 




0.1250 


0.1290 


0.1385 


0.1400 


0.1400 




0.1330 


0.1330 


0.1420 


0.1450 


0.1450 








0.1455 




0.1500 


h-t 2 


15-28 


20-28 


15-28 


20-28 


20-28 




1.88(2) 


2.55(1) 


1.88(2) 


2.55(1) 


2.55(1) 


a- l (M p ) 


1.92(11) 


2.56(21) 


2.19(9) 


2.88(14) 


2.92(20) 




2.00(6) 


2.73(18) 


1.91(6) 


2.95(14) 


2.96(11) 


a-^if^-lp-planeQ 


2.00(10) 


2.7(1) 


2.25(6) 


3.0(1) 


3.0(1) 



Table 1: Summary of the parameters of the runs analyzed in this work. 
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in eq. (|3|) and from discretization errors of 0(a), a being the lattice spacing, present 
in the operator matrix elements. The use of chiral Ward identities allows a non- 
perturbative determination of Za [f|, thus eliminating this source of error. Errors of 
order a can be reduced by the use of "improved" lattice actions @, g|. Another method 
to get rid of Za consists of extracting the decay constants of heavier pseudoscalar 
mesons by computing the ratio Rp = fp/ f n and multiplying Rp by the experimental 
value of the pion decay constant. Hopefully, by taking the ratio, some of the 0(a) 
effects are eliminated. These effects are expected to be more important for fjj than for 
f W) K since, at current values of a, the relevant parameter m c harmO is not very small. 

In the last few years, a wide set of results for fp, obtained in the quenched ap- 
proximation by using several lattice quark actions and from numerical simulations on 
different lattice volumes and at different values of the lattice spacing, have appeared in 
the literature ||. Some of these studies tried to extrapolate the decay constant to the 
continuum limit a = by using values of fp obtained at different values of the lattice 
spacing. If successful, this extrapolation would eliminate the error due to discretiza- 
tion. The extrapolation, however, is difficult because either one has to use results 
obtained at large values of a, where the dynamics is very different from that of the 
continuum limit, or the range of a is too small to observe and correct the 0(a) effects 
with a sufficient precision for the extrapolation. For these reasons, the conclusions of 
these studies are not very convincing and have changed in time. 

In this paper, we present the results of a high statistics study of fp, at two different 
values of the lattice spacing, corresponding to (3 = 6.0 and (3 = 6.2, obtained with 
the Wilson and the SW-Clover actions ]7). We have preferred to concentrate our 
computational effort on two values of (3 which have been chosen: 

a) small enough to obtain accurate results on reasonably large physical volumes; 

b) large enough to avoid the dangerous strong coupling region, which sits at around 
= 5.7. 

The main parameters of the numerical simulations are given in table |. A com- 
parison of the results from two different actions allows us to study the reduction of 
discretization errors in the improved case and to verify the validity of some recipes 
that have been proposed to correct 0(a) effects in the Wilson case [12, [l3| 0. In the 
following, we will denote these recipes with the generic name of KLM prescriptions. 

In analyzing our results we have been particularly careful to isolate the ground 
state to avoid higher-mass state contamination. We also studied the dependence of the 
final results on the extrapolation in the heavy and light quark masses. A poor control 
of these aspects in the extraction of fp can mimic spurious 0(a) effects. 

The study of the effects which can fake discretization errors, combined with the 
high statistics of the numerical runs, leads us to the conclusion that a even higher 
statistics and a larger spread of a values are required to uncover satisfactorily the 0(a) 



2 A study of the same problem, performed by using the Ward identities to determine the renormalization 
constants of the axial- vector and vector currents, Za and Zy , as a function of the heavy quark mass, mn , 
can be found in ref. Jl4|. 
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dependence of the decay constants. Further studies at smaller values of the lattice 
spacing (corresponding to = 6.4 and 6.6) with comparable (or smaller) statistical 
errors and similar (or larger) physical volumes are required to reduce this source of 
uncertainty. 

The main physical results of our study have been given in the abstract. They 
substantially agree with recent compilations made in refs. || and J15|| . Preliminary 
results of this study can be found in ref. |l6| . 



2 Details of the analysis 



In this section, we describe in detail the extraction of the decay constants from two- 
point correlation functions and the extrapolation of the results in the heavy and light 
quark masses to the physical points (including the calibration of the lattice spacing). 
We discuss systematic errors coming from the time interval of the fit, from the method 
used to derive the decay constant, from the fit used for the extrapolation in the heavy 
quark mass or in the light quark mass, and from the choice of the physical scale (from 
the string tension a, the mass of the rho mass M p or from f n ). 

2.1 Extraction of the decay constants 

Consider the correlation function 

C AP (t) = ^2(0\A (x, t)pt(0,0)|0>, (4) 

X 

where P(x) = Q(x)j$q(x). We also introduce 

Cpp(t) = £(0|P(f,i)pt (0,0)|0>. (5) 

X 

At large Euclidean times CAp{t) and Cpp{t) behave as 

Cap® = l^ e ~ MpT/2 sinh(M P (T/2 - t)) , (6) 

C PP (t) = ^ e - MpT ' 2 cosh(M P (T/2 - t)) , (7) 

where T is the temporal extension of the lattice. We extract the raw lattice value of 
fp-i fp, using the usual ratio method 



fp"* = ( coth(Mp(T/2 - t))\ #H (8) 



\Cp P (t) y ^ 1 "/ M P 

where (...) is a weighted average over a given time interval t\-t2- Mp and Zpp are 
extracted from a fit of Cpp(t) as a function of t to the expression given in eq. (0) in 
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fp = R P x = Jp Xu x f™ P » (10) 



the same interval t\— 1% in which we compute fp U . The values of t±-t2 used to extract 
all the results reported below are given, for the different simulations, in table [l|. 
The physical value of fp is then simply given by 

_ (0\A \P(p = 0)) latt , 
* P = Mp = * p ^ ' ' ^ ' 

Alternatively, we can extract the decay constant of the heavier mesons, by normal- 
izing it to (or to fx), defined as the pseudoscalar decay constant computed in the 
chiral limit 

fp^iMp) 

fla 

J TV 

where f^ u = fp U (Mp = 0). In the following, it will be useful to introduce also 
R Ps = /j?«//Jf , where f£ u = f l ^ u (M P = M K ) and /]?** is /J? tt computed at a value 
of the light quark mass corresponding to the strange quark, m s . We have then 

fPs=Rp s XfK P - (11) 

We now discuss the choice of the time intervals given in table |l] and the size of 
the systematic error due to higher mass states. In all the cases, the range was chosen 
by demanding that the contamination of the excited states in the parameters of the 
fits (ie. (§^7 coth), Mp and Zpp) is at most 20% of the statistical error. We impose 
this strict criteria to stop the systematic effects of the higher states swamping the a 
dependence of f l p U . In particular, we have verified that, due to the contamination of 
the excited states, the error on the decay constant fp, for values of the heavy and light 
quark masses close to the appropriate ones for the D s meson, i.e. without almost any 
extrapolation in run and m s , is of about 2 MeV (assuming a~ 1 (f3 = 6.0) = 2 GeV and 
a' 1 = 6.2) = 2.6 GeV) and that the relative error on R P is 5R P /R P ~ 0.005. This 
makes us confident that this error is much smaller than the other ones (the statistical 
error, the error due to the extrapolation in the quark masses and the error coming 
from the calibration of the lattice spacing). 



2.2 Extrapolation of the raw results 

In order to obtain the physical values of fp,, fp> s , etc., we have to extrapolate fp both 
in the heavy and light quark masses, to fix the value of the lattice spacing in physical 
units and to determine Za (or to use Rp to extract the decay constant). 

The best method to monitor 0(a) effects is to study Rp, computed for a fixed 
physical value of the meson mass, using the lattice spacing a calibrated by the string 
tension for which, in the quenched approximation, errors are of 0(a 2 ). By extracting 
the decay constant from Rp, we eliminate the apparent o-dependence coming from the 
variation of Za with (3. Za is in fact independent of quark masses, except for 0(a) 
corrections. Since meson masses are not exactly the same for the different simulations, 
we are forced to make some extrapolation both in the heavy and the light quark masses. 
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Method 


C60 C62 W60 W62a W62b 


linear 

Rd 3 = quadratic 
f Ds /f K linear KLM 
quadratic KLM 


1.56(3) 1.48(6) 1.11(3) 1.23(4) 1.19(5) 
1.57(4) 1.49(7) 1.13(4) 1.25(5) 1.20(5) 
1.59(3) 1.50(6) 1.48(6) 1.52(5) 1.47(6) 
1.61(4) 1.51(7) 1.51(5) 1.55(7) 1.47(6) 


linear 
Rd = quadratic 
f D /U nnear KLM 
quadratic KLM 


1.63(4) 1.58(8) 1.14(4) 1.31(6) 1.25(7) 
1.69(8) 1.73(16) 1.19(8) 1.43(11) 1.28(9) 
1.67(4) 1.60(8) 1.52(5) 1.61(7) 1.53(8) 
1.72(8) 1.75(16) 1.59(10) 1.77(13) 1.59(11) 


IdJId linear 
quadratic 


1.08(1) 1.07(2) 1.06(1) 1.07(1) 1.09(2) 
1.09(3) 1.04(4) 1.06(3) 1.06(3) 1.13(3) 



Table 2: Summary of the physical results for Rp> s = fo s / ftc an d Rd — fo/ fir obtained by 
extrapolating Rp and Rp s . We also give /d s //d- "linear" and "quadratic" refer to the fit in 
the light quark masses. 



In order to minimize the extrapolation, which will be discussed below, and at the same 
time obtain a physical quantity, we have chosen Rd s = fo s l fx-, which is obtained from 
Rp 3 []. Indeed, for all the runs listed in table |], one of the heavy quark and one of the 
light quark masses are very close to their physical values {m c harm and m s respectively) , 
and so in this case the effect of the extrapolation is negligible. Only the extrapolation 
in one of the light quark masses (to the chiral limit) is then relevant in order to get 
Rd„- It turns out, see table ||, that for Rp> a the difference between the value obtained 
with a linear or a quadratic fit to Rp s is (at most) 2 %, corresponding to an error of 
about (less) 4 MeV for fp> s . This is to be contrasted with Rp> = fo/ fir as obtained 
from a fit to In this case, differences between values obtained from a linear or a 
quadratic fit can be as large as 10 %, mostly due to a quadratic dependence of f w on 
the light quark masses. The relevance of these differences for the physical results will 
be discussed in the next section. 

As for the mass of the charmed and strange quarks, they have been determined 
by fixing the D- and if -meson masses to their physical values. The choice of the 
physical quantities used to fix the different quark masses may also affect the effective 
a-dependence of the decay constants, since the spectrum also suffers from discretization 
errors. A different possibility for the method used to fix m s will be discussed later on. 
Our results for Rpj, Rpj s and /d s //d are given in table |2| for the following cases: 
i) run C60, ii) run C62, iii) run C60 with KLM corrections, iv) run C62 with KLM 
corrections, v) run W60, vi) run W62 (a and b), vii) run W60 with KLM corrections 

3 It may also be that for this quantity the error due to the quenched approximation is smaller, although 
we have not checked this point yet. 
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(KLM- Wilson), viii) run W62 with KLM corrections (a and b). 

Since several versions of the KLM prescription exist and, in the Clover case, it must 
be implemented from 0(a 2 ) only [14], we now state the recipe used in this work to 



obtain the results in table |2|. We corrected fp U , as obtained from eq. (J8|) , for any given 
pair of values of the quark masses (m^), by multiplying it by the factor 



?klm = V(l + omi)(l + am 2 ), (12) 
in the Wilson case and by the factor [BX 



in the Clover case Q, where 

F lt 2 = l+\ 



TrC _ J ' KLM /-, o\ 

-> KLM — p p ) \ 16 ) 



(1 + amip) - (1 + amip) 



In table ^, results from linear and quadratic fits in the light quark masses are both 
given. We do not give results for fu s l jo with KLM corrections because they essentially 
cancel out for this ratio. 

We have the following observations on the results given in table a) Without 
KLM factors the results in the Wilson case are incompatible with those obtained with 
the Clover action, b) KLM- Wilson results are compatible (indistinguishable) with the 
Clover results at = 6.0 {(3 = 6.2). c) Within the statistical errors which are of the 
order of 3% (corresponding to about 6 MeV of error for fD a ), KLM- Wilson results do 
not exhibit any appreciable a-dependence; d) Even with such small statistical errors, it 
is impossible to decide whether there is a systematic shift of the Clover results between 
(3 = 6.0 and (3 = 6.2 or the difference only comes from a statistical fluctuation. This is 
due to several reasons: differences between KLM- Wilson and Clover results at f3 = 6.0 
are of the same size than those between (3 = 6.0 and (3 = 6.2 in the Clover case; in the 
Clover case, the KLM corrections to Jd s / Ik go in the wrong direction (they increase 
the difference between the results at (3 = 6.0 and 6.2). If we fit Rp linearly in the 
light quark masses, the KLM corrections increase the difference between the results 
at (3 = 6.0 and 6.2 also in the case of Ro- This is opposite to what happens with a 
quadratic fit, although within larger errors: in this case the KLM corrections reduce 
the differences between the results at (3 = 6.0 and 6.2. e) In the Wilson case, the 
variations observed between two different runs at (3 = 6.2 are about one half of the 
differences between the two Clover results (3 = 6.0 and (3 = 6.2. 

From the above considerations, it is clear that any attempt to extrapolate our 
results to a = in order to reduce discretization error would be fruitless. This is 



even more true because, as observed in ref. [11|, the results of the extrapolation are 



4 



The factor \/ AK\ K2, where K\$ are the hopping parameters corresponding to the masses mi^, is 
always included in our definition of the currents. 
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extremely sensitive to the choice of the scale, given the small range in a at disposal. We 
could have easily done a high statistics, very accurate calculation at lower values of (3, 
e.g. /3 = 5.7, to increase the range in a. We believe, however, that at such large values 
of a the behaviour of the lattice dynamics is very different from the continuum one, 
so that the inclusion of this point could fake completely the final result. This is true 
not only for our data, which have tiny statistical errors, and have been obtained on 
(approximatively) the same physical volume, but also for many other analyses which 
have been recently presented in the literature ||. With this we do not mean that in 
the KLM- Wilson and Clover cases there aren't sizeable 0(a) effects for mu around the 
charm quark mass ^ (indeed there are indications, from calculations of the dependence 
of the effective Za on run , that these effects may be large [14Q. We want to stress that, 
in spite of the very good accuracy of our data, it would be illusory to try to correct 
0(a) effects by extrapolating to a = 0. This can be done only by enlarging the range 
of a towards even smaller values, corresponding to /3 = 6.4 or even 6.6, while keeping 
the physical volume, and the statistical accuracy the same as at f3 = 6.0 and 6.2, see 
also ref. flTi 



2.3 Other sources of uncertainty 

In this subsection, we discuss some tests which we performed in order to check the 
stability of the results for the D-meson decay constants (see table |2|) and the B-meson 
decay constants (presented in table ||). These calculations consist in extracting the 
physical results with different assumptions on the parameters and the method chosen 
for the extrapolation. 

To be specific, we consider Rd s , obtained in the Clover case from a linear fit in the 
light quark mass and (then) from a fit in the heavy quark mass (at fixed m s ) of the 
form 

feV^^"^^ + -§ + ■■., (14) 

S / ,s 

and similarly for fp^/Mp. <3?* n ^, & s , are functions which are expected to depend 
logarithmically on uih but have been taken constant in the fit. For Rp> g the inclusion 
of the quadratic term of eq. (|i~4] ) in the fit is immaterial, since it amounts to a difference 
of less than about 3 MeV for fp> g . In the -B-case, the difference is of about 5 MeV 
for Jb s - Equation (|TJ) has also been used to extrapolate in mn the results of table 
||, discussed in the previous subsection. If not stated otherwise, the fits in the heavy 
quark mass always include the quadratic term. 
We have considered the following cases: 

1. Choice of the scale. Although Rp> s is a dimensionless quantity, the calibration of 
the lattice spacing can affect its value because it enters in the determination of 
the values of the quark masses at which we extrapolate Rp s - At f3 = 6.2 with the 

5 It is evident that without the KLM factors these effects are very large in the Wilson case. 
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Clover action, using a linear fit in the light quark masses, we obtain Rd s = 1.48(6) 
(scale from a), 1.48(6) (scale from M p ), 1.47(5) (scale from /„. with Za = 1.045 as 
determined non-perturbatively with the Ward identities ]18[) and Rp> s = 1-47(6) 



(scale from K*, using the lp-plane method of ref. [11], see below). A similar result 



is obtained at (3 = 6.0 where we used Za = 1-06 [19|. Thus the error due to the 



calibration of the lattice spacing is negligible for this ratio. 

2. Strange quark mass: There are different equivalent methods to fix m s , from Mk, 
from Mk*, or with the "lattice physical plane" (lp-plane) method which has 



been used in refs. [11, £(]] and that we briefly recall below. All of these methods 
should give the same value of m s , apart from discretization errors and quenching 
effects. Thus we also computed Rjj g by fixing m s with the lp-plane method for 
comparison. With this method, for the problem at hand, one defines a physical 
plane [My a, (Mpa) 2 ], where My is the mass of the vector meson. In [My a, 
(Mpo) 2 ], assuming that only linear terms are important, one looks for the point 
where My/Mp coincides with the physical value Mk* /Mk, and reads off the 
value of m s which is then used to evaluate Rd s - The results that we obtain in 
this case are indistinguishable from those given in table 

3. Extrapolation in the heavy quark mass. Among the suggestions introduced with 



the aim of correcting 0(a) effects in the Wilson case, it has been proposed [21], 
besides rescaling the quark fields according to the KLM prescription, to shift the 
mass Mp by 

M P a^ Mpa = M P a-m H + m H , (Q,q) -> M(K H ,i) (Q,q) , (15) 
where 

A K 

M(K H ) = _ 3)1/2 ^ rh H = log{Af 2 {KH)h 

m,H = e ma sinh mjj/(sinh m# + 1) . (16) 

K cr n is the critical value of the Wilson hopping parameter and N{Kh) is simply 
a different form of the KLM prescription which was used in eq.(|i~2"D. We have then 
fitted the result in Mp instead than Mp. The values of Rp> s change as follows: 
R Ds = 1.48(6) 1.44(3) (W60), 1.52(5) -> 1.51(5) (W62a) and 1.47(6) -> 



1.47(6) (W62b). Thus, the results obtained by using the recipe in eq.(15) are 
indistinguishable, within the errors from the KLM- Wilson results reported in 
table §. This is not surprising because Mp is only shifted by a small amount from 
Mp and the range of the extrapolation is very small. The effect is only slightly 
larger after the extrapolation to S-mesons, and give a small shift, corresponding 
to an uncertainty of about 5 MeV for fp s ■ We included this error, by combining 
it in quadrature with other ones, in the final evaluation of fp 3 and fp, which can 



be found in subsection 3.2 
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This concludes the discussion of several minor uncertainties in the calculation of Rd s , 
fD s /fD, Rb s and /b s //b- All of them have little influence in both the D- and 5-meson 
cases. 



3 Physical results 

On the basis of the discussion of section we are now ready to present our final results 
and errors. We will first give the results for charmed mesons, for which the extrapola- 
tion in the heavy quark mass is not a relevant source of systematic uncertainties, and 
then discuss the S-meson case. 



3.1 Results for fp s and fp 

From section ^, we learned that an extrapolation to a = it is not possible at this 
stage. Thus we believe that the best estimate of fn s is obtained from the Clover data 
at (3 = 6.2, by using Rd s in eq. ( |TT[) (from a linear fit in the light quark masses, a 
quadratic fit in 1/Mp s and without any KLM factor). As for the error, we take as 
a conservative estimate of the discretization error the difference between the results 
obtained at /3 = 6.0 and 6.2, and combine it in quadrature with the statistical one. 
This gives Rd s = 1.48(10) from which, by using f^ p = 159.8 MeV, we obtain fo s = 
237 ± 16 MeV. By using R Ds = 1.48(10) combined with f D Jf D = 1.07(4) we obtain 
fo = 221 ± 17 MeV. We believe that these results, which have also been given in the 
abstract, are our "best" results. They are in very good agreement with the compilation 
of lattice calculations presented in refs. || |I5| . Since at (3 = 6.2 the ratio Rd s is 
essentially identical in the Clover and KLM- Wilson case, the Wilson results do not 
add much information, besides giving a consistency check. 

If we extract /b from Rp, by using eq. ([h]) instead, we have to take into account 
the larger differences due to the use of the linear or the quadratic extrapolation in 
the light quark masses. We also take into account an error on /b of about 6 MeV, 
which comes from the difference in the results obtained with a linear or quadratic fit of 



fp to eq. (14). Proceeding as before, i.e. taking as discretization error the difference 
between the result obtained at f3 = 6.0 and 6.2, and combining it in quadrature with 
the statistical error, we find f D = 209 ± 16 MeV (linear) or f D = 228 ± 22 MeV 
(quadratic). Our previous result for fjj sits in the middle of these two numbers. 

To reduce the uncertainties due to the extrapolation in the light quark masses, and 
which are mainly related to the quadratic dependence of on the quark masses, we can 
compute directly Jd, by assuming some value for the calibration of a and for Za- By 
taking a -1 as determined from M p , which in the Clover case is a _1 (/3 = 6.0) = 1.92(11) 
and a -1 (/3 = 6.2) = 2.56(21), Z A (j3 = 6.0) = 1.06 and Z A (f3 = 6.2) = 1.045 |f, 
and fixing m s from Mr, we obtain (only linear fits in the light quark masses) /b(/9 = 
6.0) = 212 ± 20 MeV, f D {(3 = 6.2) = 204 ± 20 MeV, f Ds ((3 = 6.0) = 219 ± 17 MeV 
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Method 


C60 C62 W60 W62a W62b 


linear 

fB s /fK quadratic 
linear KLM 
quadratic KLM 


1.48(7) 1.28(9) 0.79(4) 0.83(4) 0.84(5) 
1.49(9) 1.27(12) 0.81(5) 0.84(5) 0.84(5) 
1.56(7) 1.33(9) 1.29(4) 1.26(6) 1.16(7) 
1.57(12) 1.33(11) 1.33(6) 1.27(8) 1.16(7) 


linear 
fs/fn quadratic 
linear KLM 
quadratic KLM 


1.53(9) 1.32(13) 0.81(5) 0.86(6) 0.86(6) 
1.57(16) 1.37(66) 0.87(9) 0.91(11) 0.86(9) 
1.61(10) 1.38(13) 1.32(6) 1.31(8) 1.19(8) 
1.65(32) 1.43(50) 1.41(11) 1.39(15) 1.19(11) 


fB a lfB linear 
quadratic 


1.10(3) 1.14(6) 1.05(2) 1.10(3) 1.12(3) 
1.13(7) 1.17(17) 1.03(6) 1.14(6) 1.20(7) 



Table 3: Summary of the physical results for Jb s I Ski /s/A obtained by extrapolating Rp 
and Rp s . We also give /b s //b- "linear" and "quadratic" refer to the fit in the light quark 
masses. 

and /d(/3 = 6.2) = 230 ± 9 MeV, which are slightly lower than the previous results, 
but perfectly consistent with them. 

The latter two determinations of fo s and fu are, however, subject to larger sys- 
tematic effects, either due to the choice of the fit or to the assumptions for the values 
of a and Za, than those extracted from Rp> s and /d s //d- We therefore prefer these as 
our best values. 

3.2 Results for fs s and /g 

In order to obtain fp 3 and fp, an extrapolation in the heavy quark mass well outside 
the range available in our numerical simulations is necessary. Discretization errors 
can affect the final results in two ways. Not only do they change the actual values of 
the decay constants, but also deform the dependence of fp on m#. Moreover, points 
obtained at the largest values of m#a become the most important, since we extrapolate 
in the direction of even larger values of ma- 
in fig. p], we show the KLM- Wilson and Clover results for fp/ fn\] /cr 1 / 2 as 
a function of the dimensionless scale a 1 / 2 /Mp. We note the remarkable agreement 
between the scaled KLM- Wilson and the Clover data (as was first observed in pHq] ). In 
the figure, we do not show the Wilson results without the KLM prescription, because 
they are inconsistent between each other (at the two values of (5) and with the Clover 
results. 

The results reported in table || were obtained by fitting Rp s (Rp) to eq. (OJ), 
including the quadratic term. We notice that all the differences observed in table || for 
D-mesons are amplified when we extrapolate to i3-mesons. Not only is this true for 
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Figure 1: Dependence of fpj ' f^[Mpj V 1 / 2 ) 1 / 2 on a 1 / 2 /Mp. For these points a linear extrapo- 
lation in the light quark masses to the chiral limit has been used. 

the results of the runs C60 and C62, where we can suspect that the differences are due 
to discretization errors, but also for the runs W62a and W62b, performed at the same 
value of (3. Thus, the same general considerations made in section ^ for Rp> s apply 
here. In order to extract our "best" values for f Bs and fp we proceed exactly as in 
subsection [H]. We obtain f Bs = 205 ± 15 ± 31 MeV= 205 ± 35 MeV, where the second 
error (31 MeV) is the "discretization" error, estimated by comparing the results from 
C60 and C62. We also obtain fp = 180 ± 32 MeV. As in the D-meson case, these 
results are in good agreement with the compilation of lattice calculations presented in 
refs. |, 0]. 

4 Conclusion 

Using several runs obtained using the Wilson and Clover actions at f3 = 6.0 and 6.2, 
from a careful analysis of all possible effects which can fake discretization errors, we 
conclude that, using the methods outlined above, we require even higher statistics, 
or a larger spread of a values, to uncover satisfactorily the 0(a) dependence of the 
decay constants. We believe that other studies of the same problem have the same 
difficulties as us in controlling and correcting discretization errors. Further studies, 
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with comparable (or smaller) statistical errors and physical volume, at smaller values 
of the lattice spacing, corresponding to (3 = 6.4 and 6.6, are required to reduce this 
source of uncertainty. The use of the action of ref. Q can be of great help in this 
respect. 

By assuming quite conservative discretization errors we found 

f Ds = 237 ± 16 MeV , f D = 221 ± 17 MeV , 
f Ba = 205 ± 35 MeV , f B = 180 ± 32 MeV , 

4^ = 1.07 ±0.04, ^ = 1.14 ± 0.08 ( 17 ) 

JD JB 



in good agreement with previous estimates 0, 15] 
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